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Abstract 

We investigate ac Josephson effect in the resonant tunneling through meso- 
scopic superconducting junctions. In the presence of microwave irradiation, 
we show that the trajectory of multiple Andreev reflections can be closed 
by emitting or absorbing photons. Consequently, photon-assisted Andreev 
states are formed and play the role of carrying supercurrent. On the Shapiro 
steps, dc component appears when the resonant level Eq is near the position 
Vl + \m\bj = Vr + (e = h = 1), where mi and ni2 are integers, uj is 

the frequency of microwave, Vl and Vr the chemical potentials of the super- 
conducting leads. Analytical result is derived in the limit A — > oo, based on 
which new features of ac Josephson effect are revealed. 
PACS numbers: 74.50.+r, 73.63.Kv, 73.23.Ad 
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Introduction-Since the achievement of atom-size superconducting quantum point con- 
tact (SQPC) by break junction technique [1-4], considerable theoretical works have been 
devoted to the transport problem in S-I-S or S-N-S structures (S=superconductor, I=tunnel 
barrier, N=normal metal). Comparison between experimental data and theoretical calcula- 
tion suggests that the coherent picture of multiple Andreev reflection (MAR) is the central 
conception in understanding the transport through such mesoscopic SQPC. Nevertheless, 
less attention is paid to the resonant tunneling (especially in ac transport) through meso- 
scopic superconducting junctions, in which the central region is consist of discrete electronic 
states. Such resonant tunneling can be achieved in nanoparticles sandwiched between su- 
perconducting membranes [5], or metallic clusters absorbed on the STM substrate [6], or a 
piece of carbon nanotube tunnel coupled to superconducting electrodes [7], or quantum dots 
fabricated in the S/2DEG hybrid structures [8], etc. 

The transport through these structures is greatly modified due to the existence of discrete 
energy levels in the central region. It is shown by Yeyati et al. [9] and Johansson et al. 
[10] that the combination of MAR processes and resonant transmission gives rise to a rich 
subharmonic gap structure (SGS). Let us consider the transport through a resonant level 
with width V coupled to S leads with superconducting gap A. In the limit T ^> A, the 
resonance is sufficient broadened such that the I-V characteristics mimic that of SQPC. In the 
regime T ~ A, pronounced SGS appears with concomitant negative differential conductance, 
which is distinguished from the nonresonant transport through SQPC. In the limit Y <C 
A, although SGS is more particular in the log-scale plot, the tunnel current as a whole 
is exponentially small. This is because a narrow resonance can not cover all the MAR 
trajectory, and it is very unlikely that multiple discrete levels are just in the position where 
MAR trajectory passes (see Fig. la). 

In this paper, we address the effect of microwave (MW) irradiation on the resonant 
tunneling through mesoscopic superconducting junctions. We show that the tunnel current 
is greatly enhanced on the Shapiro steps, when the resonant level E is near the position 
Vl + \m\io — Vr + |m 2 u; (e = h = 1 throughout), where m\ and m 2 are integers, u is the 
frequency of MW, V L and V R are the chemical potentials of left and right S lead. This can 
be attributed to the formation of photon-assisted Andreev states (PAAS), which play the 
role of carrying supercurrent. In the limit A — > oo, analytical results are derived, revealing 
new features of ac Josephson effect in the case of resonant tunneling. 

Hamiltonian and Formulation-The model Hamiltonian reads H = Y1b=l r + H cent + 
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H T , in which Hp = Ylka £ k a \ka a Pka + J2k(A el ^ a pk-\ a p-ki + H- c ) is tne Hamiltonian for 
the /3th S lead, H cent = E J2 a c \ c <? the resonant level in the central region, and H T = 



i JlVf}{t')dt' a^ ka c a + H.c.) the tunnel couplings. Vp(t) = Vg + VFgcoswt is 
the time-dependent voltage drop across the /3th barrier, where Vp is the chemical potential 
controlled by dc bias voltage and Wp cos cot the ac voltage induced by MW irradiation. For 
simplicity, we assume that Coulomb interaction and multi-level effect can be ignored in the 
central region. This assumption is somehow too ideal to achieve, but allows us to obtain 
analytical results and therefore instructive for the understanding of more complicated cases. 
Keldysh Green functions are defined in the 2x2 Nambu representation: 



/ir,o,< 
^1*2 



^ «CT(*i)i4(*=)» r,tt,< << c t(*i)M* 2 )>r' a ' <A 



V 



/ 



((^ tl )\c\(t 2 )))^< ((clWlcfoW" 
The time-dependent current flowing out of the /3th lead can be expressed as 
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Ip{t) = 2ReTra z (G r o £< + G< o 
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in which o z is the 3rd Pauli matrix, o the shorthand notation for integration over intermedi- 
ate time variable, and the self energy due to tunnel coupling between the central region 
and the /3th lead. G r , G a , and G K satisfy the integral equation: 
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E|(c) = /(c)[S^(c)-Ej&(e)], (9) 
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with being the coupling strength, rj the dephasing rate in the S lead, /(e) the Fermi 
distribution function, and Im y/z > as a convention. The remaining task is to solve these 
integral equations and evaluate the dc component of the time-dependent current. 

There are two intrinsic frequencies in the problem, u>i = 2V = 2{Vl — Vr) and u 2 = cu. 
Generally, one may perform a Fourier transform 

and derive the recursive relation for A W2 (e), as done in [11]. The calculation in this way 
relies on the numerical computing power, and analytical results are impossible. We note, 
however, that in the case of narrow resonance, dc component appears only if uj\ = Nuj 2 with 
N being an integer. This becomes clear by considering PAAS shown in Fig. lb and lc: an 
electron (a hole) incomes through the resonant level E and Andreev reflected by the right 
superconducting gap, the reflected hole (the reflected electron) exchanges mi photons with 
the MW field so that it can again passes through Eq. Then the hole (the electron) is Andreev 
reflected by the right superconducting gap as an electron (a hole), and exchanges m 2 photons 
to close the trajectory (detailed discussion for photon-assisted Andreev reflection is available 
in [12]). It is easy to see that the formation of PAAS requires 2V = (mi + m 2 )u;. For this 
reason, we shall only consider the case of V — Vn — fw, while the current deviated from 
this condition is negligible small. 

The problem is largely simplified since u> can be used as the basic frequency in the Fourier 
expansion. Define the Fourier transformation as 

A lt2 =J2^ luJtl [ ^e-^-^Me), (12) 
i J 71 

A mn (e) = i m _ n (e - nu). (13) 

The definition guarantees the nice property that if C = A o B then C m „(e) = 
J2k A m fe(e)Bfe n (e). The Fourier transformed g r ' a and E^' a,< are 

C(«)=Wl£-H (14) 
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J m - l - N (a L )E r £^i(e l )Ji- n {a L ) J m _2(a L )£^<(e+) J n _ ; (a L ) 

in which ap = ^ is the MW strength on the /3th tunnel barrier, J n (x) the nth Bessel 
function, e° = e — luo and = e — /c<j ± Vjv, Vl — V/v and Vr = are set as a convention. 
Correspondingly, the equations for G r , G a , and G K are Fourier transformed into G r ' a (e) = 
g r ' a (e) + g r ' a (e)S r ' a (e)G r ' a (e) and G<(e) = G r (e)£<(e)G a (e), or equivalently, 

G r ' a (e) = g r ' a (e) + g r ' a (e)£ r ' a (e)g r ' a (e) + g r ' a (e)£ r ' a ( e )g r ' a ( e )£ r ' a (e)g r ' a (e) + ■■■, (17) 
G<(e) = [g r (e) + g r (e)^(e)g r (e) + • • •] £<(e) [g a ( e ) + g a ( e )£ a ( e )g a (e) + •••]. (18) 

We note that finite order perturbation expansion is inadequate in the problem, because 
the formation of PAAS involves up to infinite order of tunneling processes. To re-sum up 
the series, we adopt the resonant approximation [12] 

1 1,1 



(e - hu ±E + i0+) (e - 1 2 uj ± E + i0+) ~ 5hh (e - hu ± E + i0+) 2 ' ^ 
The approximation implies that the overlap between sidebands E + liu and E q + 1 2 uj can be 
ignored if li ^ l 2 , which is justified when Fp <^ cu. Applying the approximation to Eq.(17) 
and Eq.(18), one can obtain the solution 
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It can be shown that the relations G r = (G a )^ and G< = — (G < ) t still hold in the solution. 

With these Green functions, the time-dependent current can be expressed in a summation 
over ac components 

l p (t) =2Re^ e ^|| ^-Tra z [G-(c)E|(e) + G<(e)EJ(e)] |0 } . (25) 

The current formula can be applied to the study of ac harmonics. However, we are more 
interested in the dc component I — II — —Ir, due to the experimental reasons. To produce 
analytical results, we take the limit A — > oo, due to which all single particle processes are 
forbidden and Andreev reflection is the only conducting mechanism. After some algebra and 
take the limit ?7 — > (notice that lim^o x2 ^2 = k5(x) ), one can derive 

J = 2s in0^(-) Jm- N (2a L )J m (2a R ) * [F(E+) + F(E~)} , (26) 



+ (27) 



F(e) =(e + E - ^mu?j ^ / + ^mu?j jf(a R ) + / + ^mu^j jf(a L ) 
- E + ^mu?j f($- \ mu ^j J -i( a R) + / (V ~ \ rm ^j J ~i( a L) 

in which £± = ±yj (E - \muf + |rj 2 and T m = \ \V R Jj2a R ) + T L e~^ J m _ N {2u L )\. 
Eq.(26) is the central result of this paper, which is for the dc component of Josephson 
current in the resonant tunneling through mesoscopic superconducting junctions. 

Before numerical study, we make a few remarks on this result: (1) The phase dependence 
of / is mainly determined by the prefactor sin </>, and a weak cos <fi dependence is hidden in 
E^ via Y m . For this reason, <j> is set as | in the numerical calculation. (2) For the special 
case N = 0, o>l = a R = 0, and T L = T R = T, I = 2sin^i [/ (-£ ) - / (+#o)] , 

with Eq = \J ' Eq + (r cos |) 2 , which reproduces the exact result for dc Josephson current 
through a resonant level. (3) The current formula is for the gauge choice of V L = V N and 

V R = 0. It is easy to see that the formula is invariant under the transformation L < ► R, 

V N — * —V^, N — > — N, and E — > E —V N , which corresponds to the gauge choice V L — and 
V R = —Vn- (4) E^ are the poles of G r (e), which can be interpreted as PAAS. Obviously, E^ 
contribute to the supercurrent with opposite signs. One can expect that resonant structures 
will appear near E = ^muj. (5) Bessel functions enter the current formula not in the 
square form Jn(x). It can be shown that I(aL,ot R ,(f)) = I(—aL,—a R ,(fi) for A" even and 
I{otL, ocr, 0) = I(—otL, —ocr, 4> + tt) for A" odd. 
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numerical results-Firstly, we discuss the case that MW is applied symmetrically to the 
left and right tunnel barrier, i.e., ai = cur. Fig. 2 shows the curves of I c = 7(0 = |) vs 
E — ¥k at bias voltage Vn = fw, with N = 0,1,2,3,4 (the curves are shifted by ^ L 
for demonstration). Four features are noteworthy in the plot: (1) The curve is symmetric 
(anti-symmetric) with respect to E = Vl + Vr = for N even (odd). (2) Photon-assisted 
structures appear near E — Vl — j(N — m)uj — Vr + \mio. (3) There are two types of 
resonant structures, single peak (dip) and peak-dip pair. (4) The structures grow with the 
MW strength non- monotonously. . Feature (1) is due to the relation I(Eq, 0) = I(Vn — Eo, 0) 
for N even and I(E , 0) = I(V N - E , + it) for N odd. Using I(E , -0) = -I(E , 0), one 
can obtain I(E , = |) = (—) N I(V N — E , = |). Feature (2) can be understood in terms 
of PAAS shown in Fig. lb and lc: Electron and hole are Andreev reflected back and forth 
by the superconducting gaps. When E is near the position Vj, — |(iV — m)uj — V R + \muj, 
quasiparticles may exchange N — m photons with the MW field at the left tunnel barrier 
and exchange m photons at the right, so that the trajectory is closed and bound states 
are formed. It is these PAAS that play the role of carrying supercurrent. Understanding 
of feature (3) and (4) needs quantitative analysis of each photon-assisted structure. To 
proceed, we decompose the total current into I = Yl m ^mi and expand I m near E = |mu;. 
Let Eq = ^mco + S, one can obtain I, 
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The cancelation between J+ and I m results in two types of resonant structures: C\ , — ^= = + 

V<5 2 +|r m | 
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The oscillatory nature of Bessel functions in these coefficients is responsible for feature (4). 
There is an interesting special case, N = 2, m — 1, oll — ctR, due to which C\ = C<i = 
for arbitrary MW strengths. Correspondingly, the resonant structure near E — ^ = is 
always missing for AT = 2. 
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Next, we investigate the case that MW is applied only to one of the tunnel barriers, i.e., 
ol = and an ^ 0. Fig. 3 shows the curves of I c vs E for ai = and or = 1. In contrast 
to the symmetric case, a single peak is pinned at E = Vl- The reason is as follows: since 
no MW is applied to the left barrier, Andreev tunneling through this barrier occurs only 
when E lines up with Vl, while photon-assisted processes is allowed at the right barrier 
irradiated by MW field. The inset shows I c vs MW strength an at Eq = Vl for N > 
(The case of N < can be easily deduced from N > 0). One can see in the plot that the 
peak height (including the sign) is predominated by the prefactor J m _Ar(2ai) J m {2ctR). For 
«l = 0, J m _Ar(2a£ / ) requires m = N, and the peak height is proportionate to Jjv(2af?). We 
note that there exist some regions of MW strength where I c and Vn have opposite signs, 
which can be viewed as quantum pump effect. However, this feature is dramatically different 
from the quantum pump effect in normal mesoscopic junctions. 

conclusion-To sum up, we have investigated ac Josephson effect in the resonant tunneling 
through mesoscopic superconducting junctions. We show that PAAS play an essential role 
for conducting supercurrent through a narrow resonance, in which MAR trajectory is closed 
by exchanging photons with MW field. On the Shapiro steps TV = Nuo, dc component 
appears when the resonant level E is near Vl — (N — m)uj — Vr + mu, due to the formation 
of PAAS. In the limit A — > oo, analytical result, Eq.(26), is derived for the dc component of 
Josephson current, and help to understand new features of ac Josephson effect in the case of 
resonant tunneling. We have dropped Coulomb interaction in the calculation, the derived 
results are meaningful in the following senses: (1) They are directly applicable to the systems 
where Coulomb blockade effect is negligible, i.e., U ~ T <C A (U is the charging energy). 
This is possible by using substrate with large dielectric constant to reduce U or using proper 
material as S lead to obtain large A. (2) They are instructive for more complicated cases. 
We note that the conception of PAAS is also useful for the interacting case. For U ^> T, 
the resonant level E is effectively split into two resonances E and E + U. Similar to 
Fig. lb and lc, one can draw diagrams of closed MAR trajectory through these resonances. 
Moreover, Coulomb blockade can be partially removed by bias voltage or MW irradiation, 
as long as U is comparable to A. Obviously, calculation including interaction term is much 
more difficult, and analytical results are almost impossible. Efforts along this line are still 
in progress. 
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FIGURE CAPTIONS 

Fig. 1 Schematic diagram of the resonant tunneling through mesoscopic superconducting 
junctions, (a): Without MW irradiation, MAR trajectory can not be covered by a 
narrow resonance, and the tunnel current is exponentially small, (b) and (c): In the 
presence of MW irradiation, MAR trajectory can be closed by emitting or absorbing 
photons. Two types of PAAS are formed, carrying supercurrent with opposite signs. 

Fig. 2 The curves of I c = 7(0 = § ) vs E - \ at bias voltage V N = f w (N from to 4), 
with symmetric MW strengths ai = olr. Parameters are: u — 1, V L — Y R — 0.02, 
k B T = 0.001. V L = V N and V R — are set as a convention. 
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3 The curves of I c = I(4> — f ) vs E at bias voltage Vjv = fw (iV from —4 to 4), with 
MW strengths ol = and ckr = I. Other parameters are the same as Fig.2. The inset 
shows I c at E = Vl vs the MW strength ctR. Diagrams for corresponding PAAS is 
also shown in the plot. 
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